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ABSTRACT 

The distance of the Fermi-detected blazar 7-ray emission site from the supermassive black hole is a 
matter of active debate. Here we present a method for testing if the GeV emission of powerful blazars 
is produced within the sub-pc scale broad line region (BLR) or farther out in the pc-scale molecular 
torus (MT) environment. If the GeV emission takes place within the BLR, the inverse Compton 
(IC) scattering of the BLR ultraviolet (UV) seed photons that produces the 7-rays takes place at the 
onset of the Klein-Nishina regime. This causes the electron cooling time to become practically energy 
independent and the variation of the 7-ray emission to be almost achromatic. If on the other hand 
the 7-ray emission is produced farther out in the pc-scale MT, the IC scattering of the infrared (IR) 

MT seed photons that produces the 7-rays takes place in the Thomson regime, resulting to energy- 
dependent electron cooling times, manifested as faster cooling times for higher Fermi energies. We 
demonstrate these characteristics and discuss the applicability and limitations of our method. 

Subject headings: galaxies: active — quasars: general — radiation mechanisms: non-thermal 

gamma rays: galaxies 


1. INTRODUCTION 

Blazars dominate the extragalactic 7-ray sky 
(Abdo et al. 2011a), exhibiting flares in which the 
Fermi 7-ray flux is seen to increase up to several times 
the pre-flare level (e.g. Abdo et al. 2010b). Because the 
emitting regions are unresolved, the location of the GeV 
emission remains unknown and this constitutes an open 
issue in understanding energy dissipation in blazars. 

In the leptonic model, the GeV emission of blazars is 
produced by IC scattering of photons off of the same 
relativistic electrons in the jet that give rise to syn- 
chrotron radiation (for a review see Bottcher 2007). The 
seed photons for IC scattering are synchrotron photons 
originating within the relativistic jet (synchrotron - self 
Compton, SSC scattering) (e.g. Maraschi et al. 1992) or 
photons originating external to the jet (external Comp- 
ton, EC scattering), such as UV accretion disk photons 
(Dermer et al. 1992), UV photons originating from the 
BLR (Sikora et al. 1994), and IR photons originating 
from molecular torus (MT) (e.g. Blazejowski et al. 2000). 

Because the 7-ray emission of blazars can vary sig- 
nificantly on short times, down to the instrument lim- 
its, ranging from a few minutes (in the TeV range) 
(Aharonian et al. 2009) to a few hours (in the GeV 
range) (Foschini et al. 2011), light-time travel arguments 
constrain the size of the emitting region and it has been 
argued (e.g. Tavecchio et al. 2010) that this requires the 
emission to take place at the sub-pc scale. Explain- 
ing the Fermi GeV spectral breaks seen in powerful 
blazars (Abdo et al. 2010c) as a result of pair absorp- 
tion on the He Lyman recombination continuum and 
lines, produced in the inner part of the BLR, places the 
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GeV emission well within the BLR (Poutanen & Stern 
2010). Very short distances < 10 16 cm), however, where 
accretion disk photons dominate are not favored, be- 
cause of the high pair production 7-ray opacity (e.g. 
Ghisellini & Tavecchio 2009). 

Gamma-ray flares have been associated with radio and 
optical flares, placing the GeV emission outside the BLR: 
Lahteenmaki & Valtaoja (2003), using the temporal sep- 
aration of the radio and 7-ray fluxes, argued that the 
7-ray emission site is at a distance of ~ 5 pc from the 
central engine. Connections have been also drawn be- 
tween 7-ray outbursts, optical flare polarization changes, 
and ejections of superluminal radio knots on the pc scale: 
Marscher et al. (2010) associated short 7-ray variations 
with the passage of disturbances from the radio core, 
parsecs away from the black hole. In the case of OJ 287, 
Agudo et al. (2011) found 7-ray flaring to be simulta- 
neous with mm outbursts associated with a VLBA jet 
component 14 pc downstream from the radio core. 

Here we propose a method to distinguish if the GeV 
emission of powerful blazars takes place within the sub- 
pc BLR or further out within the pc-scale MT. In §2 
we discuss our understanding of the BLR, MT, and SSC 
photon fields and where each one dominates. We show 
that cooling on the UV BLR photons results in achro- 
matic cooling times in the Fermi band, while cooling 
on the MT IR photons results to faster cooling times at 
higher 7-ray energies. In §3 we present our method, dis- 
cussing the dilution of an energy-dependent cooling time 
from light travel time effects and discuss the applicability 
of our method. Finally, in §4 we present our conclusions. 

2. COOLING OF GEV EMITTING ELECTRONS 
2.1. Characterizing the BLR and the MT 

Reverberation mapping of the BLR points to a typi- 
cal size of Rblr ~ lO 17 !^^ cm (e.g. Kaspi et al. 2007; 
Bentz et al. 2009), where 45 is the accretion disk lu- 
minosity in units of 10 45 erg s _1 . Because of this scaling, 
and assuming that L B lr = £ blrLcl , where £blr ~ 0.1 
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is the BLR covering factor, the energy density inside the 
BLR ( Ublr ~ 10 -2 ergcm -3 ) is not expected to vary 
widely from object to object (Ghisellini & Tavecchio 
2009). As was shown by Tavecchio (2008), when trans- 
formed to the jet frame, the BLR spectrum, dominated 
by H Ly a photons, is broadened and boosted in lu- 
minosity due to relativistic effects and can be well- 
approximated by a blackbody spectrum that peaks at 
Rblr = 1-5 T v Lya ( cblr « 3 x 10 _5 r), where T is the 
bulk Lorentz factor and e is the photon energy in me 2 
units. 

In the MT the dust temperature ranges from T ~ 300 
K (Landt et al. 2010) (to T ~ 1200 K (Cleary et al. 
2007; Malmrose et al. 2011), with hotter emitting ma- 
terial expected to be closer to the central engine (e.g. 
Nenkova et al. 2008). Because of the larger size of 
the MT, it has been probed by reverberation mapping 
only for Seyfert galaxies (e.g. Suganuma et al. 2006) 
which have lower luminosities and therefore smaller MT 
sizes. These studies, along with near-IR interferomet- 
ric studies (e.g. Kishimoto et al. 2011) also support an 

R oc L y 2 scaling, suggesting that U is similar for 
sources of different luminosities. Adopting the scaling 

of Ghisellini & Tavecchio (2009), Rmt = 2.5 x 10 
cm and assuming a covering factor £mt ~ 0.1 — 0.5, 
the MT photon energy density is Umt ~ 10 -4 erg 
cm -3 . We also adopt a temperature Xmt ~ 1000, corre- 
sponding to a blackbody peak at vmt ~ 6.0 x 10 13 Hz 
(«mt ~ 5 x 10~ 7 ). 


2.2. Which seed photons dominate where? 

For an isotropic photon field, the co-moving (jet frame) 
energy density U' is U' ~ (4/3)T 2 t/, while for pho- 
tons entering the emitting region from behind, U' ~ 
(3/4 )T~ 2 U (Dormer & Schlickeiser 1994). If the emis- 
sion site is located at R < Rblr , the BLR photon 
field can be considered isotropic in the galaxy frame 
and, using the BLR energy density discussed in §2.1, 
U'blr ~ 1.3 x r 2 0 erg cm -3 . Similarly, the MT pho- 
ton field is isotropic inside the BLR and its co-moving 
seed photon energy density is U' MT ~ 1.3 x 10 -2 T 2 0 
erg cm -3 . Clearly, at Rblr scales, U' BLR dominates 
over U' MT by a factor of ~ 100. If the emission site 
is located at R ~ Rmt , then the BLR photons en- 
ter the emitting region practically from behind, so that 
U' B lr ~ 7.5 x 10~ 7 T/q 2 erg cm -3 . The MT photons re- 
tain the same co-moving energy density as before and, 
in this case, the MT dominates the co-moving photon 
energy density. 

These external photon field co-moving energy densities 
need to be compared to the synchrotron one. If Rbiob is 
the size of the emitting region, then the co-moving syn- 
chrotron photon energy density is U' s ~ L s / (47rci? 2 /ob <5 4 ), 
where <5 is the usual Doppler factor and L s the syn- 
chrotron luminosity. An upper limit to Rbiob is set by 
the variability timescale t v : Rbiob ct v 8 /{ l + z), where 
z is the source redshift. We then obtain a lower limit 
U' s « L s ( 1 + z) 2 /(47rc 3 f 2 <5 6 ) or U' s ~ 6.3 x 10 _2 L Sj 4 6 (1 + 
z:) 2 t~^ d ^io er § cm_3 i where t„,id is the variability time 


in days. Setting S = T the condition U' s < U' ext requires 


/ 3L s (l + z) 2 y _ / £ M6 (l + ^) 2 \ i 

\lbTrcH 2 v U' ext ) ' \tl,i d u Lt,-i ) 


(1) 


This translates to T > 4.8 x L^ 4 8 6 (l + z) 1 ^H v 1 /f for 
the emission site at Rblr scales and V > 8.6 x L s46 (l + 

for the emission site at Rmt scales. In general, 
powerful quasars exhibit superluminal motions requiring 
T - 10 - 40 (Jorstad et al. 2005; Kharb et al. 2010) and 
therefore, the dominant source of seed photons for these 
sources should be the BLR or the MT, provided the y- 
ray emission site is at a distance not significantly larger 
than the size of the MT. A recent finding by Meyer et al. 
(2012) that the ratio of 7-ray to synchrotron luminos- 
ity increases with increasing radio core, as expected for 
EC emission (Dermer 1995; Georganopoulos et al. 2001), 
supports further the possibility that powerful blazars are 
EC emitters. 


2.3. Cooling in the BLR vs Cooling in the MT 

The critical difference between the BLR and the MT 
is the energy of the seed photons: the BLR produces 
UV photons while the MT produces IR photons. This 
difference by a factor of ~ 100 in typical photon energy 
is critical in that it affects the energy regime in which the 
GeV-emitting electron IC-cooling takes place, and thus 
the energy dependence of the electron cooling time. 

For this consideration to be relevant, IC cooling must 
dominate over synchrotron cooling. This condition is sat- 
isfied, since in powerful blazars the IC luminosity clearly 
dominates over the synchrotron luminosity, reaching in 
high states an IC luminosity higher that the synchrotron 
one by a factor of up to ~ 100 (e.g. Abdo et al. 2010a; 
Meyer et al. 2012). 

For electrons cooling in the Thomson regime (yeo << 
1, where both the electron Lorentz factor 7 and the seed 
photon energy e 0 are measured in the same frame) the 
cooling rate 7 oc y 2 . For electrons with yeo >> 1 cool- 
ing takes place in the Klein-Nishina (KN) regime with 
7 oc In 7 (Blumenthal & Gould 1970). For the broad in- 
termediate regime, 10~ 2 < yeo < 10 2 , a parametric ap- 
proximation has been suggested (Moderski et al. 2005). 
In what follows we calculate numerically the electron en- 
ergy loss rate 7, which for monoenergetic seed photons 
of number density no and dimensionless energy eo is 

7 = no [ ~T~~f( x )i e ~ e o)de, (2) 

Jo 4 eo 

where, following Jones (1968), 

f(x) = 2x log x + x + 1 — 2x 2 + 1 (4e °7 :r) (1 — x) (3) 

1 + 4e 0 yx 

with x = e/[4y 2 eo (1 — e/y)]. 

The effects of the transition between Thomson and KN 
regimes on the electron energy distribution (EED) and 
the resultant spectrum of the synchrotron and IC emis- 
sion have been studied before (e.g. Blumenthal 1971; 
Zdziarski 1989; Dermer & Atoyan 2002; Sokolov et al. 
2004; Moderski et al. 2005; Kusunose & Takahara 2005; 
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Fig. 1. — Electron cooling time in the galaxy frame as a function 
of 7 . Top panel: Blazar emission site located in the BLR. Bottom 
panel: Blazar located in the MT. Broken lines represent the various 
cooling mechanisms and the solid black line is the total cooling 
time. Plots were calculated for seed photon energies eg blr = 
3 X 10 ~ 5 and eo ,MT = 5 X 1 0 “ ' , seed photon energy densities 
Ublr = 1 0 2 erg cm -3 , Umt = 10 — 4 erg cm ' 3 and magnetic 
field energy Ub = 100 Uec- 

Georganopoulos et al. 2006; Manolakou et al. 2007; 
Sikora et al. 2009). 

In short, the transition from Thomson to KN cool- 
ing in the EED can be seen for the steady state 
77.(7) c* / Q{t)^T / T- for a power law injection <5(7) oc 
7 -p , the cooled EED hardens from 77(7) oc 7VP+1) 
in the Thomson regime to 77(7) oc 7 _ ( p_1 ) in the KN 
regime. The transition takes place gradually at 7 ~ 
Cq 1 . A second transition back to 77(7) oc y _ ( p+1 ) is 
expected at higher 7, as synchrotron losses (assumed 
to be less significant than Thomson EC losses) be- 
come larger than KN losses. The imprint of these con- 
siderations on the observed synchrotron and IC spec- 
tra has been considered for the case of the X-ray 
emission of large scale jets (Dernier & Atoyan 2002) 
and in the case of blazars (Kusunose & Takahara 2005; 
Georganopoulos et al. 2006; Sikora et al. 2009). 

Spectral signatures, however, are not unique and can 
be diluted or altered by extraneous causes (e.g. the 
accelerated electron distribution has intrinsic features 
that deviate from a power law, the synchrotron emission 
is contaminated by the big blue bump, the IC emission 
is modified by pair production absorption from syn- 


chrotron or external photons). For these reasons their 
use in understanding where the GeV emission takes 
place is problematic. 

3. LOCATING THE GEV EMISSION SITE 

There is however another observable aspect of the 
Thomson - KN transition that is free of the issues faced 
by relying on spectral signatures: the energy dependence 
of the electron cooling time can be used to evaluate the 
regime at which the electrons that produce the GeV radi- 
ation cool, and through this evaluate where the emission 
takes place. 

In the Thomson regime, the electron cooling time 
t c = 7/7 scales as 7 -1 , while in the KN regime the 
electron cooling time scales as as 7/ In 7. These two 
regimes connect smoothly forming a wide valley around 
760 ~ 1 in which the cooling time is practically energy- 
independent. The energy of the seed photons originating 
in the BLR is greater than that of photons from the MT 
by a factor of ~ 100. As a result, (see Fig. 1) the lo- 
cation of the energy-independent valley of the BLR case 
is manifested at electron energies lower by a factor of 
~ 100 than that of the MT case. In turn, the energy 
dependence of t c causes a related energy dependence in 
the observed falling time tf of the flare. The falling time 
tf of Fermi light curves can be used, therefore, to deter- 
mine whether a flare occurs within the BLR or within 
the MT. 

3.1. Flare decay in the BLR vs flare decay in the MT 

We utilize an one-zone code to demonstrate the ef- 
fect of a flare occurring within the BLR versus a flare 
occurring within the MT. The code assumes the injec- 
tion of a power law electron distribution and follows 
its evolution taking into account the KN cross section, 
SSC losses, and assuming a black body external photon 
field. The implicit numerical scheme is similar to that 
used by Chiaberge & Ghisellini (1999) and Graff et al. 
(2008). After the code reaches a steady-state, an in- 
crease in the injection rate is introduced that lasts for 
a short time before it returns to its initial value. The 
code was initialized to simulate the flare occurring within 
the BLR and the MT respective, described by the val- 
ues discussed in §2. We assumed a Compton dominance 
L ic ,max / L S ynch,max ~ 50, source size R — 3 x 10 cm, 
bulk Lorentz factor T = 10. Deep in the Thomson 
regime (MT case) the Compton dominance corresponds 
to Uec/Ub , where Uec is the external energy density 
and Ub is the magnetic energy density. However, as 
cooling transitions from the Thomson to the KN regime 
(BLR case), the ratio Uec/Ub over-estimates the ob- 
served Compton dominance by a factor of ~ 2. 

In the BLR case (Fig. 2, blue curve) IC cooling oc- 
curs at the onset on the KN regime, resulting in a near- 
flattening of tf within the Fermi band, where a change 
in energy by a factor of 100 (from 200 MeV to 20 GeV) 
results in a change in cooling time of less than an hour, 
with all the Fei'mi energies having practically indistin- 
guishable decay times (between ~ 2 — 3 hours). In the 
MT case (Fig. 2, red curve) cooling occurs mostly in 
the Thomson regime and tf is heavily energy-dependent; 
the same change in energy by a factor of 100 results in 
a change in cooling time by ~ 10 hours. This is also 
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Fig. 2. — l/e falling time in the galaxy frame versus observing 
energy (GeV) for MT seed photons (red), BLR seed photons (blue) 

apparent in the light curves shown in Fig 3. A flare 
occurring within the BLR (Fig 3, top panel) has compa- 
rable values of tf within the Fermi energy range. A flare 
occurring within the MT exhibits distinct falling times 
that decrease with observing energy (3, lower panel). 

For sufficiently bright flares with short decay times 
(comparable to the cooling times of few to several hours 
anticipated), Fermi light curves can be generated at 
multiple energies, and the energy dependence (or lack 
thereof) of tf can be used to reveal the location of the 
GeV emitting site. 

3.2. Light Travel Time Effects 

Because the GeV emitting region is not a point source, 
the observed l/e decay time t a b s of the light curve is 
a convolution of tf and the blob light-crossing time ti c . 
Any other variations (e.g. a gradual decrease of the elec- 
tron density) will have a similar effect. For BLR achro- 
matic cooling, we expect the light curves to remain achro- 
matic after considering light crossing time effects. To 
see the effects of ti c in the case of MT cooling, where, 
tf ex e" 1 ' 2 , we consider a source with exponentially de- 
caying emission coefficient, j(t) = jo exp _t 'V , and calcu- 
late the expected light curves for two energies differing 
by 10 and for a range of ti c . As can be seen in Fig 3.2, 
although t 0 b s increases with increasing ti c , the difference 
of t 0 bs between cle and Che is practically preserved for 
all tic- Moreover, because for given cle and Che the ra- 
tio oitf is also given, observations of the l/e decay times 
for cle and £he can be used to find both f; c and tf for 
each energy. 

3.3. The feasibility of the diagnostic 

Very bright flares, such as the November 2010 flare of 
3C 454.3 (Abdo et al. 2011b) have enough photon statis- 
tics to produce light curves in two Fermi energy bands. 
The practical application of our diagnostic requires that 




t (hr) 


Fig. 3. — Light curves at 200 MeV, 2 GeV, and 20 GeV for BLR 
seed photons (upper panel) and MT seed photons (lower panel), 
normalized for a steady state L=0 and maximum L=1 (in arbitrary 
units). Times are in the galaxy frame. The dashed black line is 
the l/e decay time. Times are in the galaxy frame. 



Fig. 4. — Exponentially decaying light curves of a flare in the 
MT. High and low energy luminosities (ehe / e LE = 10) are plotted 
(normalized to 1) for increasing t[ c . Solid lines represent ehe and 
dashed lines ele- Colors represent ti c at differing factors of tf, 
where tf,HE is normalized to 1. The dashed black line is the l/e 
decay time. 
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the decay times in different energy bands have statis- 
tically meaningful differences if the flare occurs within 
the MT. To evaluate if this is the case for a flare as 
bright as that of 3C 454.3, we adopted the maximum 
observed errors in the 0.1—1 GeV range and > 1 GeV 
range (Abdo et al. 2011 b) and produced simulated data 
assuming energy- dependent cooling in the form of expo- 
nentially decaying light curves. We fitted an exponen- 
tial decay function to the simulated data for each energy 
band, with the characteristic decay time as a free pa- 
rameter. The decay times we recovered were statistically 
distinguishable, suggesting that our diagnostic is appli- 
cable, at least for the brightest flares observed. 

3.4. Constraints from upper limits of the decay time 
difference 

For flares from which no statistically significant differ- 
ence in the decay time can be established, constraints 
can still be imposed on the location of the Fermi de- 
tected emission by considering A t max , the maximum de- 
cay time difference allowed by the data, between a high 
energy che and a low energy cle- The requirement that 
A t m ax > tf (c-le) — tffeHE) can be casted as upper limit 
for the location R of the Fermi emission 

2 1 1 ^ 2 

R < 2(JTLMTr 2 At max 

3-\/3 7 rm e c 2 (l + ^) 1 / 2 e^(e ^ /2 - e^ 2 ) 

. ( 4 ) 

where Lmt is the MT luminosity. In terms of quantities 
used in §2.1 and for cle, e le corresponding to 100 MeV 


and 1 GeV, R < 2.3 x 10 18 Fio(Ai TOaXi /i.LMT, 45 /(l + 
z) 1 / 2 ) 1 / 2 cm. The Lorentz factor can be estimated from 
VLBI studies (e.g. Jorstad et al. 2005) and the MT lumi- 
nosity can be taken to be a fraction £ ~ 0.1 — 0.5 of the 
accretion disk luminosity or in some cases can be directly 
measured (e.g. Landt et al. 2010; Malmrose et al. 2011). 

4. CONCLUSIONS 

We presented a diagnostic test for determining if the 
GeV emission in powerfull blazars originates from within 
or outside the BLR. Our method utilizes the fact that if 
electrons cool via IC scattering on BLR photons, cooling 
occurs at the onset of the KN regime and the resulting 
cooling times (and thus light curves) should be achro- 
matic. If cooling occurs via IC scattering of less energetic 
MT photons, the electron cooling times, and thus light- 
curves, should exhibit significant energy-dependence and 
decrease as photon energy increases. The energy depen- 
dence of the cooling time difference perseveres in the 
presence of energy independent timescales such as the 
light crossing time. The method is applicable for the 
brightest Fermi 7 -ray flares, and can provide upper lim- 
its on the flare location even when the decay times in 
different energies are statistically indistinguishable. 
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